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Abstract. A polaron model proposed as a possible mechanism for nonlinear 
conductance [Galperin M, Ratner M A, and Nitzan A 2005 Nana Lett. 5 125-30] 
is revisited with focus on the differences between the weak and strong molecule-lead 
coupling cases. Within the one-molecular level model we present an approximate 
expression for the electronic Green function corresponding to inelastic transport case, 
which in the appropriate limits reduces to expressions presented previously for the 
isolated molecule and for molecular junction coupled to a slow vibration (static limit). 
The relevance of considerations based on the isolated molecule limit to understanding 
properties of molecular junctions is discussed. 
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1. Introduction 



Much of the interest in molecular conduction junctions stems from their functional 
properties as possible components in molecular electronic devices. In particular, non 
linear response behaviours such as bistability and negative differential resistance (NDR) 
have attracted much attention. Here we revisit a model for such phenomena that was 
previously advanced [1] and later criticized [H [Sj H] in order to eludidate and clarify some 
of its mathematical characteristics. 

The simplest molecular conduction junction comprises two metallic electrodes 
connected by a single molecule. The simplest theoretical model for such a junction 
is a molecule represented by one electronic level (the molecular affinity or ionization 
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level) with one vibrational mode connecting free-electron metals. When the molecular 
electronic level is outside the range between the lead Fermi levels and its distance from 
these levels is large compared to the strength of the molecule-lead electronic coupling, 
the transport occurs by tunneling through the molecular energy barrier. This is the so- 
called Landauer-Imry limit. When the injection gap (distance between the Fermi level 
and the affinity or ionization levels) becomes small, the barrier decreases, and there is an 
opportunity for stabilizing excess charge on the molecule by polarization of its electronic 
and/or nuclear environment, leading to the formation of polaron-type trapped charge. 
We have previously described the consequences of this polarization on such phenomena 
as hysteresis, switching and negative differential resistance in molecular junctions. [Ij 

When the electronic coupling between the molecule and leads vanishes, one deals 
with polaron formation on an isolated molecule, for which an exact solution is available. 
We discuss here the two limiting cases: polaron formation on an isolated molecule, 
and the transport problem in the limit where nuclear dynamics is slow relative to all 
electronic timescales. Invoking the second case as one of the possible mechanisms of 
hysteresis, switching, and negative differential resistance in molecular j unctions [1] was 
criticized by Alexandrov and Bratkovsky, in several papers. [21 |3l H] These authors claim 
that the conclusions of Ref. P contradict a previously published "exact solution" [5l [6] 
that shows no multistability is possible for molecular models comprising nondegenerate 
and two-fold degenerate electronic levels. They suggest that multistability found in 
Ref. [Ij is "an artifact of the mean-field approximation that neglects Fermi-Dirac 
statistics of electrons" {hq = no), and "leads to a spurious self- interaction of a single 
polaron with itself and a resulting non-existent nonlinearity" . 

As was pointed out previously, [7] the weakness of this criticism stems from using, 
in Ref. [5], the isolated quantum dot limit to discuss molecular junctions. In contrast, 
we have argued [7] that the approximtion of Ref. [1] is valid in the limit T ^ ujq, where 
F is the inverse lifetime of excess carrier on the bridge and ujq - the frequency of the 
relevant nuclear motion. Here we present this argument in a rigorous mathematical 
form. We describe a general approach to this problem, which is capable reproducing 
the result of Ref. [5] in the isolated molecule limit and our previous result, Ref. [1], in 
the static limit of a junction (cuq/F -C 1), where uq is the oscillator frequency and F, 
the spectral density associated with the molecule-lead coupling, measures the strength 
of this coupling. This validates the polaronic approach of Ref. [1] in this limit. 

2. General consideration 

One way to bridge between the limits of zero and strong molecule-lead coupling is 
the nonequilibrium linked cluster expansion (NLCE) proposed in Ref. [8]. For our 
purposes a first order LCE[10] (clusters of second order in electron-phonon coupling 
Mjjl is adequate. Indeed, this level of consideration provides exact results in both 
isolated molecule and static limits, while providing an approximate expression for the 

t We use the term "phonon" for any relevant molecular or environmental vibration. 
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general case. The main idea of the NLCE is the same as in the usual LCE - one 
expands a Green function (GF) perturbatively in terms of the interaction part of the 
Hamiltonian (in our case - the electron-phonon interaction) up to some finite order, 
and equates the expansion in clusters to an expression in terms of cumulants. [9] This 
provides approximate resummation of the whole series. [TOj The NLCE considers this 
expansion on the Keldysh contour p| 



T, r 



n=0 



T, T 



Ln=l 



(1) 



whence, up to first order [n = 1) 

Woir, t') = Go(r, r') W^ir, r') = Go(r, r')Fi(r, r') (2) 

Projections of ([1]) on the real time axis are obtained using Langreth rules [TTl [12], in 
particular 



ln=l 



(3) 



G>'<(t,t') 



(4) 

W>'<{t,t')=G>'<{t,t')F,>'^{t,t') (5) 

In steady state (which we consider below) projections depend on time difference t — t' 
only. 



3. Model 



As in Ref. [T] we consider a single (nondegenerate) electron level Eq coupled to one 
vibration uq and to two leads L and R represented by reservoirs of free electrons, each 
in its own equilibrium. The vibration is represented by a free oscillator at thermal 
equilibrium. The Hamiltonian of the system is (here and below e = 1, m = 1, and 
h = l) 

H = eod^d + Y. {^kclck + Vkd^Ck + V^cld) 

k£{L,R} 

+ LUod^d + M{d + d))d)d (6) 

where d {d)) and Ck {c\) are annihilation (creation) operators for electrons on the 
molecule and in the contacts respectively, while d (d^) are annihilation (creation) 
operators of a vibrational quantum. The first and second terms in ([6]) represent electrons 
on the bridge and in the contacts, respectively and the third and fourth terms describe 
molecule-leads coupling. The fifth term describes the free vibration, while the last is 
the linear electron-phonon coupling. For future reference we also define the operator of 
molecular level population 

fiQ = d)d (7) 

and its quantum and statistical average 

/+00 (j TP 
^G<{E) = -iG<{t = Q) (8) 
-oo 27r 
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where G< is the electron lesser GF.jHl [12] 



4. Mathematical evaluation of transport properties 



The non-equilibrium Green function technique provides a convenient framework for 
evaluating the desired transport properties. To obtain the steady-state current under 
given bias conditions 



P r+oo 

Ik = t dtTi [S< (-t) G>{t) - S> (-t) G<it)] 

ll J —CO 



(9) 



{K = L, R) one needs to evaluate the molecular electronic Green function in the presence 
of the moleule-lead and electron-phonon couplings. In what follows we derive this 
expression within the low-order NLCE described in Section O 

The free phonon GFs (retarded, advanced, lesser and greater) are 
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where Nn 



1] ^ is the thermal equilibrium vibration population. 



In the absence of electron-phonon coupling, M = 0, electron GFs in the wide band 
approximation (where the spectral densities F;^ = lirY^ki^K |Vfcp5(-E — ^k) are energy 
independent) are 
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(14) 
(15) 

(16) 

(17) 



Tk [K = L, R) are the electron escape rates from the molecule due to coupling to left and 
right leads, F = F^ + F/j and fxiE) = [e^^~^''^^'^ + 1]~^ is the Fermi-Dirac distribution 
in the contact K {fi^ is chemical potential). In approximations made in f|T6l) and f|T7j) 
we have used uq ^ (Fl/l(£:o) + ^nfRi^o)) /F. Note that these approximations are used 
for convenience only and do not influence the generality of the considerations below. 
They become exact either in the case of molecule weakly coupled to contacts or when 
molecular level is far (compared to F) from the contacts' chemical potentials. 
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Figure 1. Lowest (A/^) order contributions to electron GF {Wi): (a) Hartree and 
(b) Fock (Born) terms. The wavy line represents free phonon GF. The straight line 
represents electron GF. 



The lowest order in the electron-phonon couphng (M^) contribution to the 
electronic GF is given by 

W^i(r, r') = j dn! dT2 G'o(r, n) r^) Go(r2, r') (18) 

J c J c 

Sph(ri,r2) = 5{Ti,T2)M^n^ jdT^Do{n,T^) +iM^Do{ti,T2)Go{ti,T2) (19) 

where self-energy (SE) is a sum of two contributions: the first and second terms in 
Eg. (1191) are respectively the Hartree and Born terms shown in Fig.[Tl The importance of 
including the Hartree term when considering systems without translational periodicity 
(e.g. molecular junctions) was emphasized in a number of papers. [131 Ull E] 

The lesser and greater projections of ( ITSj) onto the real time axis (here and below 
we assume steady-state situation) are obtained from the Langreth rules [TTl [T2] as 

r' + OO 



dti 



dU 



Gl{t-t,)^X{t,-t,)Gl{t, 



(20) 



+ ^0' (t - ti) Sp^(ti - t2) ^0(^2) + ^^^(t - ti) Sp/j(ti - ^2) Gq ' (^2) 
Projecting (fT9|) and using Eqs. (fT0l) - (fT71) one gets 
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It should be emphasized that the uq term that enters the Hartree contribution in 
Eqs. ffTI?]) . fl^Tl) and fl2^ is an exact result; unrelated to the convenient approximation 
made in Eqs. (fT6|) and (ITTl) above (that leads to the explicit appearance of the no terms 
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in Eqs. f l23p and (12^ ). It is this term which will provide the population dependent shift 
of the electronic level in the static limit (see below). 

Our aim is to get an expression for the retarded electron GF 

G^\t) = e{t)[G>{t)-G<{t)] (25) 

using ©-([S]). In order to do so we have to calculate Wi'^ which is given by Eq. (!20l) . It 
is convenient to consider separately the first term and the sum of the second and third 
terms on the right-hand-side in fl20l) 



wc^{t) = wc^{t) + w^'^(t) 
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utilizing ([I4])-([I7D and ([21])- (121 then leads to 
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(26) 
(27) 

(28) 
(29) 



(30) 



The last term in curly brackets on the right-hand-side in (l30ll comes from the Hartree 
term. The expression for (t) is obtained from fl29l) by interchanging A'^o and A'^o + 1 
and replacing no by no — 1. W^(t) is obtained from ( l30l) by replacing no by no — 1 
only in the prefactor that multiplies the curly brackets on the right-hand-side. These 
general approximate (first order LCE) expressions for >V^2^ are the central result of this 
consideration. 
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5. Two physical limits 

In [1] we have discussed a mean field approach to describe non-linear response of 
molecular junctions characterised by strong molecule-lead coupling as well as slow 
vibrations strongly coupled to the electronic subsystem. As noted in the introduction 
this approach was criticised in Refs. [21 |3l H] as incompatible with observations made in 
the isolated molecule. To elucidate the issue we consider next these two specific limits: 
the isolated molecule (F 0) and static limit {ojq/T —>■ 0). 



The isolated molecule 

In the limit T ^ Eqs. ([29]) and §^ yield 
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and the corresponding expressions for Wi2(t) 
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Substituting ([3lD-(lMD into ([26]) and using Eqs. ([3]) and ([5]), one gets from (|25D 
G'it) 
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where 
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Eq. (l35!) is the standard expression for the retarded Green function in the isolated 
molecule case, obtained following a small polaron (Lang-Firsov or canonical) 
transformation. |9] In particular, it is identical to Eq.(30) of Ref. [5] for the case of 
a nondegenerate level (i.e. d = 1 there). Note that approximations f[T6l) and f[T7l) 
become exact in this limit and, furthermore, the first order LCE provides the exact 
result in this limit. As was pointed out by Alexandrov and Bratkovsky[2l [S] H] the 
electronic level shift. A, is independent of level population for the isolated molecule, 
and no multistability is possible in this case. 
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The static limit 

The ciJo/r — > limit reflects either a slow vibration or a strong molecule-lead coupling. 
For molecules chemisorbed on metal and semiconductor surfaces F is often of order 
0.1 — 1 eV, so this limit is expected to be relevant for the relatively slow molecular 
motions associated with molecular configuration changes. To describe the behaviour of 
our model system in this case we expand the exponentials and the fractions in Eqs. fl29l) 
and (l30i) in powers of oJo/T, disregarding terms of order higher than 1 and keeping in 
mind that due to the e~^'*'/^ prefactor ojq/V ~ uj^t holds. This implies 

which leads to 

Wi<(t) = (38) 
W<{t) = moe-o*-ri*l/2 2^^^t (39) 

and corresponding expressions for VVf2(^) 

Wf(t) = (40) 
y^>{t) = -^[1 -no]e-^^°*-ri*l/^^^^^t (41) 

Substituting ( l38l) -( l4Ti) into ( l26i) and using the result in (j3]) and ([5]), one gets from ( l25l) 
G'{t) = _,0(t)e-*(-o-2noA)t-rt/2 ^42) 

Again we note that the factor no that enters this expression does not result from 
approximations ( JT6l) and (|T71) . Rather, it arises from the exact expression for the 
Hartree term, the first term on the right-hand-side in Eq. ffTOj) . Note also, that the 
approximation used in Eqs. ( |T6l) and ( |T7|) could in principle be relaxed. This would make 
the mathematical evaluation more difficult (unless the molecular level is far, compared 
to F, from the leads' chemical potentials, when this approximation becomes exact) but 
would not influence the estimates of W{^2^ in terms of coq/T. 

Note that technically the static limit corresponds to disregarding all diagrams (in 
all orders of electron-phonon interaction) except the Hartree term (see Fig. [1^) and 
terms of similar character (only diagrams with boson lines terminated in a closed loop), 
since these are the only diagrams transmitting zero frequency. In the static limit this 
is not a mean-field approximation but an exact result. Detailed discussion on the issue 
can be found in Ref. [T4j . 

To conclude, in the static limit (which is the limit considered in Ref. [1]) the 
electronic level shift, 2noA, does depend on level population in the way presented in our 
polaron model. [1] In what follows we briefiy reiterate the implications of this observation 
on the conduction properties of molecular junctions with strong coupling between the 
electronic and nuclear subsystem. [1] 
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6. Non-linear conduction in static limit 



Here we discuss briefly the consequences of tlie reorganization energy dependence on tlie 
average electronic population in the molecule, as presented in Eq. fj42|) . on the junction 
transport properties. Since we consider steady-state transport, i.e. all GFs and SEs 
depend on time difference only, we can go to the energy domain. The Fourier transform 
of Eq.(ll2]) is 

G'-{E) = [E-e,{n^)+iV/2] (43) 

where eo(^o) = — 2noA is the population dependent energy of the molecular level. 
Using the Keldysh equation 

G>^<{E) = G''{E) T>^<{E) G\E) (44) 
in expression (IHl) for the level population leads to 

^ r+^dE hiE)TL + fR{E)Tn , . 

i-oo 27T [E - BoirioW + [r/2]2 ^ ' 

This is the central result of Ref. yy (see Eq.(13) there). The non-linear character of 
Eq. ( I45l) with respect to no leads to possibility of multistability, and results in non-linear 
character of the junction transport. In particular, the zero-temperature case allows 
analytical evaluation of the integral. We find that Eq. (l45l) is equivalent to the following 
pair of equations (see Eq.(20) of Ref. [1]) 



= ^ arctan [x + ^) + ^ arctan (x - ^) + i 

4A 2A 

where V is source-drain voltage. System of equations (146|) defines points of intersection 
of an arctan function with a straight line, which for some set of parameters may have 
multiple solutions. Detailed discussion of the consequences of this multistability for 
transport can be found in Ref. [1] 



7. Conclusion 



In this paper we have presented solid theoretical foundations for the polaron model of 
non-linear response of molecular junctions, which was proviously introduced using mean 
field arguments. We have used the non-equilibrium linked cluster expansion to second 
order, and focused on the limit of the isolated molecular polaron on one hand, and the 
polaron formation in a functioning molecular transport junction (that is, with finite 
coupling to the electronic states in the leads) on the other. Proper examination shows 
that the former case indeed requires integral charge on the molecule (this is self evident, 
since there is no source or drain for the electrons). The functioning junction can have a 
non-integer average population of electrons on the molecule, and is maintained at steady 
state by the actual current flow through the molecule. 

This formal analysis demonstrates the validity of the polaron model as originally 
suggested, and shows clearly an example of a new molecular regime of functioning 
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transport junctions, characterized by strong molecule-lead coupling and slow molecular 
vibrations strongly coupled to the electronic population on the molecule, where the 
junction effect on its environment can be described by its non-integral electronic 
population. Furthermore it shows that in this case, due to the phonon polarization, the 
electronic level energy becomes dependent on this population. This is not a "spurious 
self- interaction" (as suggested in Refs. [2l[3lll]), but rather describes the interaction of 
a tunneling electron with its predecessor (s) via the phonon polarization cloud created 
by the electronic transient density of the molecule. 

Finally, while we believe the mathematical issues concerning the model advanced 
in Ref. [1] has now been clarified, it should be pointed out that actual observations of 
multistability and NDR in molecular junctions can arise from other mechanisms. In 
particular, to account for such observations in the Coulomb blockade regime we would 
probably need to go beyond the simple model considered here. 
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